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ABSTRACT. We present a new quantization scheme for 2D gravity coupled to an SU (2) 
principal chiral field and a dilaton; this model represents a slightly simplified version of 
stationary axisymmetric quantum gravity. The analysis makes use of the separation of vari- 
ables found in our previous work |Q] and is based on a two-time Hamiltonian approach. The 
quantum constraints are shown to reduce to a pair of compatible first order equations, with 
the dilaton playing the role of a "clock-field" . Exact solutions of the Wheeler-DeWitt equa- 
tion are constructed via the integral formula for solutions of the Knizhnik-Zamolodchikov 
equations. 



Introduction. At present only few models of quantum gravity exist for which the 
Wheeler-DeWitt equation can be solved exactly. The known examples include pure gravity 
and supergravity Q in three dimensions, as well as certain mini-superspace models such 
^Supported by Alexander von Humboldt Foundation. 

1 



as static spherically symmetric gravity or locally supersymmetric Bianchi-type models 
[^]. All of these models describe only finitely many physical degrees of freedom and are 
thus bound to miss essential features of quantum gravity. It is therefore clearly desirable 
to find models of quantum gravity with infinitely many physical degrees of freedom, which 
can still be treated exactly. In this paper we present such a model describing 2D gravity 
coupled to a principal chiral SU (2) model and a dilaton, and demonstrate that it can be 
quantized exactly by use of methods developed in the context of 2D integrable systems 
1^, ^. The construction is based on our previous work Q, where we have proposed a 
novel canonical formulation of the Ernst equation based on the complete separation of 
the equations of motion in the isomonodromic sector of the theory and the use of the 
logarithmic derivative of the related ^'-function with respect to the spectral parameter as 
the fundamental canonical variable. Our model can be regarded as a version of the "midi- 
superspace" approximation of quantum gravity, but differs from dimensionally reduced 
gravity in that we have replaced the non-compact configuration space 1)/C/(1) by 

the compact group manifold SU{2). Apart from this technical modification (the reasons 
for which will be briefly explained at the end of this paper) the model captures all the 
essential features of dimensionally reduced quantum gravity, and is a non-trivial example 
of a completely integrable model of quantum gravity with infinitely many propagating 
degrees of freedom. A further intriguing aspect of the present work is the possible relevance 
of quantum groups suggested by the link with the Knizhnik-Zamolodchikov (KZ) equations. 

Classical model. For definiteness we consider only Euclidean worldsheets with com- 
plex local coordinates (^, ^) (the extension of our results to Lorentzian worldsheets is com- 
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pletely straightforward). The worldsheet metric is 



ds'^ = h{C,^)dCd^ h = e^^ (1) 

and the Lagrangian is taken to be 

C = p[R + tTig^g-^g^g-^)) (2) 

Here /?(^, ^) G R is the dilaton field originating from dimensional reduction as explained 
e.g. in P], and the chiral field g{C,0 ^ G = SU{2) describes the matter content of the 
model (subscripts stand for partial derivatives). The conformal factor in (jl]) represents a 
further, but dependent, degree of freedom and appears in (|^) via the Gaussian curvature of 
the worldsheet; in the conformal gauge this term becomes R = (log h)^^/h. The equation 
of motion for p{^,^) following from @ reads 

Pee- = (3) 

and possesses the general solution /)(^,C) = /(^) + f{0- Assuming / ^ const we can 
perform an analytic local change of coordinates such that /(^) = S,/2i (global aspects of 
this change of variables were discussed in , but will not concern us here) . In the sequel 
we shall thus always assume 

/9(e,o = ime (4) 

and also set x{S^, ^) := Re^. We will take some liberty in referring to p and x as the "time" 
and "space" coordinates, respectively. The gauge choice (^) is analogous to the light cone 
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gauge in string theory whereby one identifies the longitudinal target space coordinate 
with the world sheet time coordinate. The remaining (and compatible) equations of motion 
read 

2% = ^—^tT{g^9~'f 2k^ = t_ltT{g^g~^f (5) 



and 

((e - 09^9-')^ + ((C - 09^9-% = 0. (6) 

As is evident from these equations, the 2D gravitational fields and the matter fields are 
coupled through the dilaton p; furthermore it can be shown that the obvious solution 
p = const of (^) would imply 9^ = 9^ = = = 0, i.e. the trivial solution for the matter 
fields, in which case the gravitational sector would become purely topological (strictly 
speaking, this argument relies on the positive definiteness of the Cartan Killing metric 
on the group G). Therefore the model possesses no non-trivial flat space limit since the 
matter fields act as sources for 2D gravity and thus distort the two-dimensional background 
geometry. 



Equation (|6|) is the compatibility condition of the following linear system [10| 



7 1 + 7 9^ _ ggg-^ 7 1 - 7 5^ _ ffgg ^ 

9? ^-^1-7^7 1-7 dC ^-^1 + 7^7 I + 7 ^ ' 



where ^ = ^(CiCi7) is a two-by-two matrix, from which the solution of (P) can be recon- 
structed, and 7 is a "variable spectral parameter" 



2 r c + e 



li^,t,w) = —{w-^ + \/{w-0{w-0\ (8) 
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with w C a constant of integration which may be regarded as the "hidden" spectral 
parameter. Notice that the poles in the 7 plane appearing on the LHS of are produced 
by differentiation of 7 with respect to ^ and ^, respectively. 

According to the approach developed in the fundamental quantities are the loga- 
rithmic derivatives of ^ with respect to the spectral parameter 7: 



*^ = ^(7)^ A(7)^E^^ (9) 



where for notational clarity we have suppressed the dependence on the coordinates (^,0- 
The two-by-two matrices Aj{^,^) are the residues at 7j = 'y{^,^;wj) in the complex 7- 
plane, with constants wj. We here consider only the isomonodromic sector, for which the 
sum (|9|) is finite; in general, however, the sum may contain an infinite or even continuous 
number of terms (if the sum is replaced by an integral) . 

The compatibility conditions of (^) and yield the following system of differential 
equations for the matrices Aj Q: 



A 



2 ^ [Ak,Aj] ^ 2 ^ [Ak,A^ 



e-e\^.(i-7.)(i-7,) c-e,^.(i + 7fc)(i+7,) ^ ' 



As emphasized in [|[] these equations are automatically compatible unlike the original linear 
system (^). The problem of solving the nonlinear partial differential equation (P) is thus 
reduced to solving a system of ordinary matrix differential equations. Once the solutions 
Aj of ( p!o|) are known, the chiral field g{£,,^) can be recovered from ^(7) by specializing the 
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spectral parameter 7 to the values 7 = ±1, viz. 



(e - Og^g-' = 2A(7)|^^i (e - Og^g-' = 2A{j)\^^_,; (11) 

again we have not displayed the dependence on the coordinates in these expressions. The 
conformal factor h{^,^) may then be computed by integration of (^). In Q we have shown 



that it is essentially the r- function associated with (7) in the sense of |11]. 

All results so far are actually valid for arbitrary SL{2, C) (and, in fact, GL{n, C)) valued 
matrices g{^,^), but we will consider only G = SU{2) in the remainder of this paper. The 
solutions of the system ( |lO|) corresponding to g ^ SU (2) and /i G R are characterized by 
the following 

Lemma 1 Let {Aj,tj;j = 1,...,A^} be invariant with respect to the involution Aj — > 
-yl] , 7j -7j (i.e. if 7^ = -7^ then Aj e su{2); if -yj = -7^ , k ^ j then 
Aj = —AI.). Then the constants of integration in ^ and ( [71^ may be chosen in such 
a way that g £ SU (2) and /i G R. 

The conditions of the Theorem 1 in particular imply that ^(1)^ = ^(—1). In the sequel 
we restrict attention to the case 7^ = — 7j (i.e. wj G R), so that Aj G su{2) for all j. 
Regularity of ^'(7) at 7 = 00 together with (^ implies a further constraint, viz. 



N 

5: A, = (12) 
3=1 



This condition is related to the asymptotical flatness of the corresponding classical solu- 
tions. All constraints are consistent with the equations of motion (|To| ) and form a closed 



algebra with respect to the Poisson bracket (^) below. 

It was also explained in that the equations (|l^) can be interpreted as a two-time 
hamiltonian system with respect to "times" ^ and ^. The relevant Poisson brackets are 
given by 

{^(7) ® ^(/x)} = [r(7 - /i), A(7) ® / + / A{fi)] (13) 

where the classical r-matrix r(j) is equal to n/7 with IT the permutation operator in 
X C^. The dynamics in the ^ and ^-directions is governed by the complex Hamiltonians 

= ^^tr (a\1)) and H^^'> = -J—ti (^^(-1)) (14) 

respectively, which depend explicitly on (x, p) through the prefactor and the spectral 
values 7j via (^). Observe that we have (ff^^-*)^ = H^^^ by Theorem 1. As expected from 
the decoupling of variables in (|l^), the flows generated by H^^'^ and and H^^^ commute, 
i.e. {H^^\ H^^^} = 0, and this ensures that the equations of motion can be consistently 
integrated. Remarkably, the dimension of the system has been effectively reduced by trad- 
ing the "space" variable x and the "time" variable p for two "time" variables ^ and ^. 
We can thus neglect the dependence on x and p for all practical purposes and regard the 
spectral parameter currents ^(7) at a fixed but arbitrarily chosen base point {xq,pq) as 
the fundamental canonical variables. 

In view of (|l^ ) we can reexpress (||) in the form 

2k^ - H^^^ = , 2k^ - = (15) 
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Below, we will have to interpret these equations as canonical constraints a la Dirac rather 
than merely as equations determining the conformal factor; upon quantization, they will 
yield the Wheeler-DeWitt equation and the difFeomorphism constraint. 



Quantization. To quantize the model, we replace the Poisson brackets (13) by com- 
mutators in the usual fashion: 



[A(7) f A{fi)] = ih[r{-/ - fi) , ^(7) ® I + I ® A{n)] 



(16) 



The matrix elements of ^(7) thus become operators acting on a Hilbert space to be specified 
below; note that on the LHS of ([l6|), we have a commutator of operators in Hilbert space 
whereas on the RHS we have a commutator of ordinary matrices. This means in particular 
that the expansion (^) is no longer valid as an operator statement, but must be reinterpreted 
as a property of the states on which A^j) acts. The expressions for the Hamiltonians (14) 
remain unchanged. We write 



/ 



^(7) 



ih 



Hi) 
2/(7) 



2e(7) 

-Hi) 



[17) 



As is well known the Hamiltonians ( |T4[ ) can be diagonalized by means of the Bethe ansatz 



method [12|. For this purpose, one considers a lowest weight state \p) labeled by some 



analytic function p = ^(7) subject to ^(7) = — p(— 7). This state is assumed to obey 
^(7)|p) = P{'y)\p) aiid f{'~f)\p) = 0- The Bethe states are then 



\p;vi, ...,vm) ■■= e{vi)...e{vM)\p) 



(18) 
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The parameters vi,...,vm rnust satisfy certain constraints called Bethe equations (M is 
the number of excitations) in order for the Hamiltonians (|l^ to act diagonally on them; 
however, we will not need these constraint equations below. The off-diagonal operators 6(7) 
and 7(7) play the role of creation and annihilation operators, respectively (interchanging 
the role of 6(7) and /(7) turns the representation "upside down"). The classical expansion 
(ph corresponds to the special choice 



N 

p{i) = E (19) 

^1 ^ - 



By fixing the number of poles in this fashion, we have restricted the full Hilbert space to 
its soliton sector, which can be represented as a direct product H^^^ = . . . (g) Hn- 
Only on this subspace can we expand the operator ^(7) in analogy with (pi), with 



ifi 

Ai = — 
^ 2 



hj 


2e, 


2/, 





(20) 



where hj , Cj, fj are the standard Chevalley generators of SU{2) obeying the commutation 
relations 

[hj, ej] = 2ej [hj, fj] = -2fj [ej, fj] = hj 



as a consequence of (pHI). By well known arguments, unitarity then requires each Tij to be 
a representation space of SU{2) with highest weight rrij E Z. We may therefore view the 
quantized soliton sector as a system of spins located at the points {x^p) = (wj,0) 
where the classical solutions of the system (^, (^ have singularities on the worldsheet 



|13|. Since the full Hilbert space is the direct sum of its N soliton sectors, its structure 



resembles that of a Fock space with N playing the role of a particle number operator (not 
to be confused with the number M of excitations in([l8|)). 

The central task is now to solve the quantum analog of the constraints ( [l5|) , which are 
nothing but linear combinations of the Wheeler-DeWitt equation and the diffeomorphism 
constraint for our model (corresponding to translations in p and x, respectively), and read 

(2k^ - $ = (2k^ - $ = 0. (21) 

where is the full quantum state (so-called "wave function of the universe" ) . Accordingly, 
we extend the phase space by the gravitational degrees of freedom 2k^ , 2k^ and the 
variables ^, ^. To proceed we must first construct an operator representation for these new 
phase space variables. From the canonical brackets given in Q we deduce the commutators 

[2k^,C] = [2k^,^]=ih , [2%,e] = [2%,^] =0 



We take 



d d 



9? ^ 

The main advantage of this choice is that by representing (^, ^) as multiplication operators 
we salvage their interpretation as coordinates (otherwise the spectral parameter 7 would 
not remain a function but become a non-local differential operator and thus very awkward 
to deal with). It is then obvious that the two equations (|2l| ) are mutually compatible 
for the same reason that their classical counterparts (|5|) are. Recall that the worldsheet 
coordinates {^,^) appear explicitly only because we have adopted the special gauge (^) 
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identifying the dilaton field with one of the coordinates. In other words, this choice of 
gauge makes the quantum state $ time-dependent through the identification of time with 
the "clock field" p. We note that this long suspected mechanism for the emergence of time 



from the "timeless" Wheeler DeWitt equation here comes almost for free (see e.g. [14] 
for a review and further references). In a covariant treatment the gauge choice ^ would 
have to be undone, and the full quantum state would be a functional of p rather than a 
function of the worldsheet coordinates. The fact that through this choice of gauge the 
field variables 2k^ and 2A;^, which originated from the conformal factor, become canonically 
conjugate to the worldsheet coordinates, is also in accord with the interpretation of 2A;g 
and 2k^ as longitudinal target space momenta suggested in ^, and can be construed as 
evidence for a "stringy" interpretation of the theory, with our treatment being analogous 
to the quantization of the string in the lightcone gauge. 
The equations (|2^) now take the form 



in^ = H^i)^ = h(0^ (22) 



where in the A^-soliton sector <I> is an 'H^^-'-valued function of {£,,(,)■ The explicit solutions 



of (I24) may be obtained by exploiting a surprising link between (22) and the KZ equations 



1 15 1 determining the A^-point correlation functions of the WZNW model, which read 



9$ 1 ^ Qjk 



-V^^-l (23) 

dlj K Ij - Ik 

with an T-^^^^-valued function l>, where = \hj ^ hk + fj (8> Cfe + ej ^ fk {j k) is a 



linear operator acting non-trivially only in Hj and Tik- In fig] (see also ||17|, [l^, |19|) it was 
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shown that the general solution of ( |2^ ) can be obtained in terms of the following integral 
representation over Bethe eigenstates 



$= / W{'y,v) \p]Vi,...,VM)dvi...dvM 
I A 



(24) 



where 



i<k<j<N l<s<r<M 



Vr — V, 



M N 
r=l j=l 



and A is a family of cycles in u-space which do not intersect the hyperplane where W 
becomes singular; furthermore, W should be single-valued on A. For the precise definition 



of A and the proof of ( plj ) see |16, 17, 18, 2C]. It is important that the parameters v,. in 



(24) are not subject to the Bethe equations. In it was shown that the state (|2^ ) is 



annihilated by the constraint (12) iff 



N 



rrii 



(25) 



Using (g) and (^) we can now verify that the KZ equations ( |23D imply that 



$=n 



(26) 



satisfies (|2^), provided the state <I> has vanishing "total spin" in accordance with (|T2|) and 
we make the identification 



K = —1/ih. 



(27) 
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Outlook. To conclude we would like to briefly expose the extra complications that one 
faces when attempting to quantize dimensionally reduced gravity, whose distinct feature 
is the non- compactness of the configuration space inhabited by the propagating degrees 
of freedom. As we already mentioned, the relevant equations of motion can be formally 
obtained from our equations (^), @ by replacing the compact SU{2) manifold by the 
non-compact coset space SU{1, 1)/U{1), after which @ becomes just the Ernst equation. 
Due to the division by the maximal compact subgroup, which is already necessary at the 
classical level in order to prevent the energy from being unbounded below, we now encounter 
new canonical constraints necessary to ensure the decoupling of unphysical negative norm 
states. In the case at hand, it turns out that, in fact, one must impose infinitely many such 
constraints corresponding to the maximal compact subgroup of the Geroch group SU{1, 1). 
Moreover, the representation spaces for each "spin" will be infinite-dimensional (unlike the 
spaces Hj above) because of the non-compactness of the group SU{1, 1). 
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